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Abstract 

We present a new method for the experimental determination of the top quark mass that is 
based upon the mean distance of travel of b hadrons in top quark events. The dominant systematic 
uncertainties of this method are not correlated with those of other methods, but a large number of 
events is required to achieve small statistical uncertainty. Large tt event samples are expected from 
Run II of the Fermilab Tevatron and the CERN Large Hadron Collider (LHC). We show that by 
the end of Run II, a single experiment at the Tevatron could achieve a top quark mass uncertainty 
of ~ 5 GeV/c^ by this method alone. At the CERN LHC, this method could be comparable to all 
others methods, which are expected to achieve an uncertainty of ~1.5 GeV/c^ per experiment. 
This new method would provide a useful cross-check to other methods, and could be combined 
with them to obtain a substantially reduced overall uncertainty. 

PACS numbers: 12.15.Ff,13.25.Hw,14.65.Ha,14.80.Bn 
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I. INTRODUCTION 

There are important reasons for measuring the properties of top quarks Q|. If there 
are new strong interactions at the TeV scale, the top quark often plays a central role in 
the corresponding theories. The large top quark mass is important in extensions of the 
standard model, particularly when mass dependent couplings are considered. New particles 
are experimentally constrained to be heavier than other fermions but could be lighter than 
the top quark. Careful studies of top quark decay products could therefore be fruitful. 
Indeed the large top quark mass and its approximately unit value Yukawa coupling to the 
Higgs particle are often cited as a hint that the top quark has a special place in the greater 
scheme of things Q]. 

The mass of the top quark, rrit, is an important Standard Model parameter that enters 
quadratically into higher order corrections to the W mass, M^y, along with the Standard 
Model Higgs mass, Mh, which enters logarithmically. Precise values of both and rrit can 
thus be used to constrain the Standard Model expectation for Mh ■ Prior to observation 
of the Higgs, this information could be used to guide experimental searches. After the Higgs 
has been observed, it can be used as an interesting consistency check. If the observed Higgs 
mass were to be inconsistent with the prediction obtained from precision measurements of 
rrit and M^y, this would indicate that nature is not governed by a Standard Model Higgs 
sector. For these reasons it is essential to measure and to the best possible precision. 

In this article we present a new method for the measurement of the top quark mass 
that is based upon the observation that the mass is correlated in a relatively strong and 
unambiguous way with the mean distance of travel of the b hadrons formed from the 
b quarks produced in top quark decays 5]. The dominant systematic uncertainties of this 
method are associated with those factors that influence the mean b hadron decay length, 
such as b quark fragmentation, b hadron lifetimes, and our ability to understand and to 
accurately model the momenta of the top quarks and their decay products. It will be seen 
that jet energy scales, the dominant systematic uncertainties in other methods, result in a 
relatively small contribution to the uncertainty in the mean b hadron decay length. As a 
result, this technique can provide an independent measurement that can be used to cross- 
check, or be combined with, results from other methods to yield a more robust, accurate, 
and precise final value for rrit. 
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The top quark was first observed nearly a decade ago in collisions of protons with anti- 
protons at a center of mass energy of jjs = 1.8 TeV during Run I of the Fermilab Tevatron 
by the CDF and DO experiments 0, l^- Using the relatively small Run I ti event samples, 
and combining the results of both experiments, the mass of the top quark was determined 
to be Q: 

mt = 178.0 ± 2.7 (stat) ± 3.3 (sys) GeY/c (1) 

This result was the fruit of a very substantial effort by large teams of physicists on both 
experiments that benefited from an innovative improvement in methodology P]. After Run 
I, there were substantial upgrades to the Fermilab accelerator complex, including an increase 
of beam energies to ^/s = 1.96 TeV, and to the CDF and DO experiments. Tevatron Run 
II began in 2001 and is scheduled to last until slightly beyond the startup of the CERN 
Large Hadron Collider (LHC). In this period, Tevatron experiments could collect top quark 
samples that are almost two orders of magnitude larger than those of Run I. These will 
enable the CDF and DO experiments to significantly improve their measurements of rrit. Each 
experiment currently expects to achieve an overall uncertainty of roughly 2-3 GeV/c^ jlol |. 

The CERN LHC will supersede the Tevatron as the world's highest energy accelerator 
when it begins operation later this decade. The LHC will collide protons with protons at 
a center of mass energy of ^/s = 14 TeV. The 7-fold increase in energy relative to the 
Tevatron means that the dominant ti production mechanism will change from quark-anti- 
quark annihilation to gluon fusion. As a result, the tt production cross-section rises from a 



calculated value of 6.7 pb 
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,Q at the Tevatron, to 833 pb lll| at the LHC. The LHC will 
also operate at higher luminosity than the Tevatron. The net effect is that the LHC will 
produce ~8M ti events per year at its initial operating luminosity of 10'^^ cm~^ sec~^. 



This number will rise by an order of magnitude for high luminosity operation. Availability of 
such large top quark samples will mean that all measurement techniques will rapidly become 
limited by systematic uncertainties, which are expected to be ~ 1.5 GeV/c^. This will also 
be true for the extraction of the top quark mass from the mean b hadron decay length. 

In this article we present estimates for the potential top quark mass resolution of this new 
method at both the Tevatron and the LHC. In Section|H]we describe in detail the aspects of 
the method that are common to its application at both accelerators. In particular, we discuss 
the important sources of uncertainty and how they are estimated. The approaches taken for 
the Tevatron and LHC are not identical, mainly due to the fact that Tevatron measurements 





FIG. 1: Feynman diagrams showing the leading order tt production mechanisms. 



will be statistically limited, even for the most optimistic projections of integrated luminosity. 
On the other hand, the large ti event rates at the LHC allow one to quickly achieve statistical 
uncertainties below 1 GeV/c^, even in the most pessimistic operating scenarios. This allows 
us to define event selection criteria that minimize uncertainties. Descriptions of the event 
selection and several other aspects of our studies that are unique to a given accelerator 
environment are postponed until sections IIII Al and IIIIBl where we also report the results 
of our studies for the Tevatron and LHC, respectively. It will be seen that variations in 
the estimated uncertainties in the two accelerator environments can be traced, among other 
things, to a difference in the dominant ti production mechanism. 

Our basic approach is to use Monte Carlo programs IJ, |l5| currently in use by Teva- 
tron and LHC experimental groups to generate tt events at the appropriate beam energies. 
We do not perform full simulations of any of the Tevatron or LHC experiments. Rather, 
we include the effects of imperfect pattern recognition and detector resolution by smearing 
Monte Carlo generated decay lengths to obtain a reasonable match to distributions seen 
in collider data and full detector simulations. In fact, we find that detector modeling has 
relatively little effect on the estimated uncertainties in rrit for this method. As described in 
section IIII A[ for the Tevatron study we consider backgrounds such as Wbb and misidenti- 
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fication of light quark jets as b jets ("mistags") in + jets and multi-jet events. For the 
LHC, as described in section IIIIBl we choose an event selection that reduces backgrounds 
to negligible levels. In addition, a jet veto is used to suppress events containing significant 
QCD radiation, thereby minimizing the uncertainties associated with Monte Carlo model- 
ing of these processes. In section IIVI we discuss ways in which improvements are likely or 
possible. Finally, we summarize our results and present our conclusions in section IVl 

The results presented below, while representative of what one could reasonably expect 
to observe, include values for mean b decay lengths in ti events that are not intended to 
exactly predict those that one would extract experimentally. For example, the kinematic 
criteria used in event selection, as well as detector specific effects, will cause variations in 
the observed mean b hadron decay length for a given top quark mass. Nevertheless, we have 
estimated uncertainties carefully and we have verified that they are robust to variations in 
the modeling of detector effects. Our projections for the uncertainties on the top quark mass 
are thus expected to be accurate predictions of what can be expected at the Tevatron and 
LHC in coming years. 



II. THE METHOD 



The Fermilab Tevatron currently collides protons with anti-protons at y/s = 1.96 TeV. 
At this energy, the dominant ti production mechanism is qq annihilation. Valence quarks 
carry ~15% of the proton momentum on average, corresponding to a constituent collision 
energy of \^ ~ 300 GeV. The ti system is therefore produced near threshold. The ti pair 
can have transverse momentum as a result of initial state radiation. This is typically below 
100 GeV, corresponding to a relativistic boost that is very near to unity: 7 — 1 < 0.04. 
In general, there is little energy available to provide the individual top quarks with large 
transverse momenta. In spite of a 7-fold energy increase to ^/s = 14 TeV, ti events at the 
LHC will be similar to those at the Tevatron in many ways. This is mainly due to the 
fact that the dominant production mechanism at the LHC is gluon fusion {gg ti). As a 
result, ti pairs are once again produced near mass threshold ~ 350-400 GeV), albeit 
with a more substantial fraction appearing at higher V^. Fig. |2l shows distributions of \^ 
underlying ti production at the Tevatron and LHC as obtained with the PYTHIA Monte 



Carlo 



li. 
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FIG. 2: PYTHIA Monte Carlo comparison of the distributions of constituent center of mass energies 
(\/^) underlying tt production at the Tevatron and LHC for mt = 175 GeV/c^. 

In the Standard Model, the top quark decay proceeds almost uniquely as t — > Wb. In 
the rest frame of the top quark, the W and b daughters have equal and opposite momenta 
of magnitude: 

P = ^ V(l - ii^w - O/mlY - A{Mwm,/m1Y (2) 

Thus, p/c 0.4 mt for large mt, and the relativistic boost of the b quark is substantial, 
75 ~ 0.4 X {iTLt/miy). As a consequence, the h quark momentum is potentially a sensitive 
gauge of the top quark mass. One could use the correlation of the h jet energy to the top 
quark mass directly. The drawback of this approach, however, is that it depends upon the b 
jet energy measurement which suffers from a jet energy scale uncertainty like other methods. 

As an alternative, consider the mean b hadron decay length, (L), which is also correlated 
with the b momentum. For a b hadron of momentum p, mass m;,, and proper life time To 
one obtains: 
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(L) = CToPj = T, 



P 



(3) 
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This expression provides the dependence on the energy since E = y {pc)"^ — {jribC^Y ~ PC. 
The key point is that (L) can be found directly without measuring E. It is simply the average 
measured distance from the primary interaction vertex to the h hadron decay vertex. It thus 
depends upon charged particle track reconstruction as opposed to jet energy reconstruction. 
Tracks are formed by fitting "hits" to helical trajectories. The hits are measurement points 
on charged particle trajectories that are typically obtained with a precision of 5-50 /im in 
semiconductor tracking detectors, such as silicon micro-strips or pixels. 

In practice, we use the average of the transverse decay length, L^y = L \ sm6\, where 9 
is the angle of the b hadron flight path with respect to the beam axis. This is necessitated 
by the fact that the net longitudinal momentum of the tt pair is not known in hadron 
collisions. The partons themselves have broad momentum distributions within the proton 
or anti-proton. The sum over the transverse momenta of all objects in the event must 
however be zero. In terms of the transverse energy Erp = E sin 6, and mass of the b hadron, 
we have: 



Equations |21 and |3] provide a simple basis for understanding the correlation between the 
mean b hadron transverse decay length and the top quark mass. There are, however, a 
large number of additional factors that influence this correlation, and therefore need to be 
understood and taken into account. These factors can be categorized as follows: 

• Factors that affect the momentum of the top quark: 

— Initial state QCD radiation from colliding partons. 

— Final state QCD radiation from top quarks. 

— Interference of final and initial state radiation. 

— Parton density functions of the beam particles. 

— Spin correlations. 

• Factors that affect the b hadron momentum and decay length: 




(4) 
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— Fraction of b quark momentum carried by the b hadron, (i.e. fragmentation). 

— Radiation from the b quark. 

— Relative proportions of b hadron species. 

— Lifetimes of the b hadron species. 

— Modehng b hadron decays. 

— b hadron identification and vertexing efficiency ("6 tagging"). 

— Dependence of b tagging on b jet momentum. 

• Background events mistakenly identified as tt. 

• Accidental b tagging of light quark jets ("mistags"). 

• Additional factors: 

— Kinematic selection of events. 

— Tracking efficiency, purity, and precision. 

— Multiple interactions in a single beam-beam crossing. 

A. Basic procedure for determination of mt and its uncertainties. 

We now give a broad overview of the basic procedures that were used in our studies, 
followed by more detailed presentations of some of the key steps. As mentioned above, the 
specific implementations of this method at the Tevatron and LHC differ in some respects. 
Details that are only pertinent to a specific accelerator environment are deferred to sections 
IIII Al and IllIBl where we describe the results of our Tevatron and LHC studies, respectively. 

1. Overview 

The procedure we follow for our studies is relatively straightforward. We create "default" 
tt samples using the PYTHIA Monte Carlo. We do not perform a full detector simulation of the 
events. We apply event selection criteria directly to the Monte Carlo generated transverse 
momenta of electrons, muons, and jets in the events. A missing transverse energy, {^t ) 
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requirement is applied to the transverse component of the vector sum of the 4-vectors of all 
neutrinos in the event. 

For a given choice of m^, we generate a large number of events. Prior to a determination 
of Lxy for a h hadron decay, we require there to be an adequate number of daughter particles 
with significant impact parameters relative to the primary vertex that can be used to form 
a secondary decay vertex. The impact parameters of these particles are smeared according 
to relevant experimental resolutions. The particles must pass typical selection criteria. The 
event selection criteria and experimental resolutions, which are typical of Tevatron and 
LHC experimental studies, are detailed in Sections IIII Al and IIII Bl If it is determined that 
the vertex could potentially be observed experimentally, then the associated L^y value is 
obtained by simply smearing the Monte Carlo generated transverse decay length of the h 
hadron. 

This procedure is repeated for different choices of m^. We also generate background 
event samples that are selected and processed in a similar manner. We combine signal and 
background events in the proportions expected for a given experiment and choice of event 
selection criteria. For each m^, the combined signal plus background sample is used to obtain 
a mean value of the transverse decay length, (L^y). These values are plotted as a function 
of rrit and they are fit to a polynomial in rrit, as seen in Fig. 01 For a given choice of beam 
energy, event selection, and b tagging algorithm, this polynomial represents our expectation 
for the correlation of rrit with (L^y) in tt events. We will refer to the polynomial as the mass 
estimator . 

To estimate statistical uncertainties, we form an ensemble of event samples in which the 
number of b hadrons from signal and background events are each randomly selected from 
Poisson distributions for which mean values correspond to the expectations for a specific 
experiment and integrated luminosity. The L^j^ values for a given sample in the ensemble 
are generated using the shape of the L^y distribution obtained for a given choice of signal 
plus background, as described above. For each sample, {L^^y) is calculated and the mass 
estimator specific to the accelerator environment under consideration is used to obtain a 
corresponding value for m^. The resulting distribution of values for the full ensemble is 
not Gaussian, in general, since the relationship of to (L^y) is not linear over large mass 
ranges as seen in Fig. El In these cases, which are relevant for small statistics experiments, 
one obtains asymmetrical statistical uncertainties. For high statistics experiments, only 
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FIG. 3: An example of the correlation of L^y with mt for the default Monte Carlo and a particular 
choice of event selection at the Tevatron. Backgrounds, as described in the text, are also included 
and represent 25% of all entries. The correlation is well fit by a third degree polynomial over the 
full mass range as shown, and it is linear over small mass ranges. Each point represents ~50,000 
b jets in Monte Carlo ti events. 

small mass ranges need to be considered and the resulting ensemble distribution for mt is 
Gaussian. In these cases the standard deviation is taken to be the statistical uncertainty. 
The uncertainties themselves are obtained by constructing Neymann freq uentist confidence 
intervals following the prescription described in §32.3.2.1 of reference (lfil |. 

To take into account uncertainties associated with the modeling of the underlying physics, 
(e.g. initial and final state radiation, and parton distribution functions), we vary the relevant 
Monte Carlo control parameters to produce alternative event samples from which we can 
again extract values of {L^y). The original mass estimator is again used to extract values 
of m^, which differ from our default values. For any particular variation, the difference 
with respect to the default, Am^, is our estimate of the uncertainty associated with this 
particular aspect of the event model. Monte Carlo parameter variations that we use are 
consistent with the current guidance provided by the authors of the Monte Carlos, or they 
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conform to current conventions used in measurements of top quark properties at the Tevatron 
or studies performed by detector collaborations in preparation for LHC operation. 

B. Lxy distributions for h, c and light quark jets and their associated uncertainties. 

We now briefly discuss how L^y is determined experimentally and the procedure used to 
obtain a reasonable approximation to L^y using Monte Carlo generated information for h 
and c hadrons. We also discuss how we treat mistagged light quark jets. This is followed 
by a discussion of the dependence of h tagging on the Et of the h jet. We then present how 
we estimate the systematic errors on {L^y) associated with the modeling of h fragmentation 
and the uncertainty in the average h hadron lifetime. We make the reasonable assumption 
that properties of h jets as measured in Z -^hb events at LEP and SLD, and high pj- fejets 
from direct production of hb at the Tevatron, apply also to h jets from top quark decays. 
This can be justified by noting that the final stage of jet hadronization is a non-perturbative 
QCD process at a scale of order Aqcd (see for instance §4.1 of jl?!). 

1. Modeling Lxy in b, c and light quark jets 

The experimentally observed L^y distribution for b and c jets is the result of a complicated 
process of particle tracking, track selection criteria, and vertex finding. While each element 
of this process contributes to the final shape of the L^y distribution to some degree, the 
final outcome can be described by two main effects. Namely, a convolution of the true 
Lxy distribution with a Gaussian resolution function, and a skewing of the distribution that 
results from the use of tracks with significant impact parameters relative to the primary 
vertex. In this section we describe in more detail how these effects arise and how we treat 
them with a Monte Carlo simulation. 

The transverse decay length is the measured distance between the point at which the 
primary beam particles collide, (the primary vertex), and the point at which the b hadron 
decays, (the secondary vertex), as projected into a plane perpendicular to the beam axis. In 
an actual experiment, a vertex is calculated as the intersection point of two or more charged 
particle tracks which have helical trajectories in the magnetic field of the tracking detector. 
Fig. m presents a schematic view of the tracking environment for a b jet from top quark 
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FIG. 4: Schematic representation of tracks in the vicinity of the beam-beam interaction for an 
event with a high energy b jet. 

decay. 

The primary vertex is reconstructed using all tracks that can be consistently associated 
with a specific, small region along the beam line. In practice, the transverse position of 
the primary vertex is sometimes taken from the beam line itself, as measured with many 
thousands of tracks, from thousands of successively recorded events, under a given set of 
operating conditions. The position of the b hadron decay is obtained by forming a vertex from 
"displaced tracks" , (i.e. tracks that have significant impact parameters relative to the beam 
axis). Typically a track impact parameter, d, is said to be significant when Sd = d/ad> x, 
where aa is the uncertainty on d, and x is typically 2.5 or 3. The tracks are also required to 
be inside a cone centered on the axis of a jet, which is defined as the line joining the primary 
vertex to the centroid of a cluster of calorimeter towers with energies above some threshold. 
For high energy b jets, a cone radius R = v^Arp~+~A0^ = 0.4 is used, where An and A0 are 
the opening angles of the cone in pseudorapidity and azimuth, respectively [l8|. 

The sources of uncertainty in the measurement of L^y include the uncertainties in the 
positions of the vertices due to uncertainties in the reconstructed helical track parameters, 
which, in turn, are affected by hit position uncertainties and detector misalignments. Other 
sources of mis-measurement are the inclusion of fake tracks, and the association of tracks 
from other vertices. Furthermore, in ti events, the b hadrons are not fully reconstructed 



12 



and so it is often the case that the tracks used to locate the secondary decay vertex are a 
mixture of tracks from the b hadron decay and sequential c hadron decay. The sequential 
c hadron has substantial momentum and in fact travels on average 1-2 mm beyond the b 
hadron decay point. Use of such vertices is necessary to obtain high b tagging efficiency 
but requires one to loosen some vertexing criteria, such as the maximum allowed value 
associated with the fit to a single vertex, but also increases the average uncertainty in L^y. 
For the case of several displaced tracks from a single decay vertex, the uncertainty on the 
transverse position of the vertex is typically 100-200 /im. This increases to ~1 mm when 
the tracks used in the vertex come from both the b and sequential c hadron decays. 

In our studies, we use smeared Monte Carlo generated information to reproduce the main 
experimental aspects of the L^y distributions for b and c hadrons. For b and c jets we use 
the 4- vectors of daughter particles with transverse momenta above 0.5 GeV/c to calculate 
impact parameters with respect to the beam line. We take into account the curvature of 
the tracks in the magnetic field in this calculation. We then associate an uncertainty with 
the impact parameter of each track using parameterizations that vary with the transverse 
momentum of the track as: 



In this expression a represents the asymptotic impact parameter resolution for arbitrarily 
high track momenta. The second term, which depends on pr, takes into account multiple 
scattering due to interactions with the detector material. The parameter f3 thus depends 
on both the amount and location of material in the tracker relative to where track hits are 
measured. The parameter a includes the uncertainty associated with the primary vertex. 
This can be quite small when the track multiplicity is large, as is often true for tt events. If 
instead one uses the position of the beam line, the relevant uncertainty is then the transverse 
size of the beam itself, which is ~ 30 /xm at the Tevatron and will be ~ 15 fim. at the LHC. 
In our studies we use different values for these parameters for the Tevatron and LHC. In 
all cases, we require that the tracks have an impact parameter significance Sd > 2.5. We 
only consider those b and c hadrons which have at least 2 significantly displaced tracks and 
we calculate an associated L^y by shifting the true transverse decay length by a random 
value obtained from a Gaussian distribution centered at the origin with standard deviation 
ai — 1 mm. This value is consistent with Monte Carlo studies of b tagging with a full 




(5) 
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simulation of the CDF detector at the Tevatron The Gaussian smearing makes it 

possible for there to be negative decay lengths, (corresponding to the unphysical case of a 
secondary decay appearing behind the primary interaction with respect to the jet direction), 
while the requirement of at least two significantly displaced tracks skews the L^y distribution 
toward positive values. The smearing and skewing just described are the most salient features 
of the Lxy distribution as it is manifested experimentally [l^ . 

We generated large numbers of events for a range of top quark masses. We also gener- 
ated large samples of background events. For each h and c jet with at least two significantly 
displaced tracks, we determine L^y as described above. We then calculated {L^y) for com- 
binations of top quark mass samples plus background samples and fit the results to obtain 
a parameterization of {L^y) versus m^. For background events we included the L^y values 
for mistagged light quark jets. Note that we obtain our normalizations for the signal and 
badcgro^d jets contribufng to the fl,>a. distribution direetly fro,„ recent experimental 
results scaled to the particular integrated luminosities and tt production rates relevant 
to our studies. 

In light quark jets, such as those from W decay in tt events or in backgrounds, it is 
possible for significantly displaced tracks to occur and lead to a false h tag. The proba- 
bility that a light quark jet will be mistagged is typically less than 1%, and is gene rally 
well parameterized by collider di-jet data as a function of jet Et and pseudorapidity |l9|. 
Mistags are due to tracks with significant impact parameters arising from several sources. 
For the most part, they are tracks that originate at the primary interaction, but appear 
to be significantly displaced because they are on the extremes of the Gaussian resolution 
distribution. True displaced tracks from particles with significant lifetimes (e.g. i^^, A°, and 
photon conversions) are also sometimes incorrectly used in h tagging when the long-lived 
particle is not properly identified. The latter occurs when at least one of the daughter tracks 
is not properly reconstructed. Finally, there are fake tracks that result from pattern recog- 
nition errors such as the association of hits from two or more particles in the reconstruction 
of a single track. Fake tracks and tracks from long-lived particles can result in reconstructed 
vertex locations distributed more or less uniformly in space out to or beyond the radius 
of the beam pipe or first measurement layer. We model these using a flat distribution at 
positive Lxy. 

Our model of the L^y distribution for mistagged jets is shown in Fig. El The distribution 
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FIG. 5: The probability distribution versus L^y for light quark jets with at least two significantly 
displaced tracks. 

is consistent with measurements and full simulation of the CDF detector fl^. It is well 
approximated by the convolution of a Gaussian centered at L^y ~ 2 mm, with symmetrically 
rising and falling exponentials in the negative and positive L^y regions, respectively, and a 
flat distribution at positive Lr^y. As mentioned above, the distribution is skewed into the 
positive L^y region as a result of the displacement significance requirement, S^. > 2.5, applied 
to the tracks. 

As a check of the sensitivity of our results to our modeling of detector effects, we repeated 
the process described above in the absence of Gaussian smearing of L^y values, and ignoring 
mistagged jets. We found that the parameterization of {L^y) versus rrit displayed a shift 
in the average transverse displacements. The slope, however, differed negligibly from our 
previous result. The scale shift means that the two cases yield different mass predictions 
for a given value of {L^y). The invariance of the slope, however, means that the estimated 
systematic uncertainties are the same in both cases. The reason that the slope does not 
change appreciably stems from the fact that the L^y distribution for any given top quark mass 
is approximately an exponential, while the uncertainties associated with the determination 
of Lj.y are Gaussian. The mass information is contained in the slope of the distribution 
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which is not changed by convolution with a Gaussian. The visibihty of the exponential, 
however, is diminished as the width of the Gaussian increases relative to the mean of the 
exponential. In the case of top quark decays, the mean decay length is roughly 5 mm while 
the Lxy resolution is typically 1 mm or less. Thus, our studies do not predict {L^y) exactly 
for a given m^, but do provide robust predictions for uncertainties on m^. 

2. Tagging efficiency versus transverse energy. 

The dependence of the h tagging efficiency on jet transverse energy, Et, is influenced by 
many factors including the choice of the h tagging algorithm and the overall performance 
of the tracking system. For simplicity, we assume a dependence on jet Et similar to that 
measured with 6 jet enriched data samples at the Tevatron Thus, we take the efficiency 
to rise linearly in the range 15 < Et < 40 GeV from a minimum of 25% to a maximum 
of 40%. It then remains approximately constant for higher values of Et- The average 
efficiency will affect the overall statistical uncertainty that can be achieved on nit for a fixed 
integrated luminosity, while the dependence on Et affects {L^y) and can therefore introduce 
a systematic uncertainty in mj. 

Studies indicate that the shape of this distribution in data is well matched by fully 
simulated Monte Carlo samples [l9|. However, the fact that the efficiency is not constant 
in the low Et region can lead to an uncertainty in nit as a result of the jet energy scale 
uncertainty. This follows from the fact that event selection criteria usually involve thresholds 
for jets in the low Et region in order to have good signal efficiency. The jet energy scale 
uncertainty thus represents a threshold uncertainty that translates into an uncertainty in 
the energy spectrum of h jets contributing to the L^y distribution. For instance, we find 
that a variation in the minimum Et threshold by an amount equal to a jet energy scale 
uncertainty of 10%, leads to a 1.09% systematic uncertainty in mt in our Tevatron study. 

Note, however, that the h tag efficiency need not be parameterized in terms of the jet 
Et- As an alternative, the total pt of the tracks involved in the tag, or associated with the 
h jet, could be used j20|. The uncertainty is then shifted from the category of calorimetry 
to that of tracking, analogous to the basic premise of the correlation of mt with {L^y). In 
this manner, the overall uncertainty due to the variation of tagging efficiency with the h jet 
Et could be reduced substantially. 
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3. Fragmentation, b hadron species, and their lifetimes 



The b quarks from top quark decays pick up light quarks from the vacuum to create b 
mesons and baryons. The b hadron carries a fraction Xf, of the initial b quark momentum, 
while the remainder is dispersed into the momenta of the lighter mesons created in this 
fragmentation process. Together with the b hadron decay products, these tracks form the 
b jet, which is generally well contained within a cone of radius R = 0.4 in tt events. The 
fragmentation function, namely the distribution of Xb values for many b jets, has been 



studied extensively. The most precise measurements o: 



of Z ^ bb events at the LEP and SLC accelerators 
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b fragmentation are from the study 



22M 



. For our studies, we use 



the single most precise measurement which comes from the OPAL experiment |2ll|: 



Xb = 0.7193 ± 0.0016 (stat) (syst). (6) 

Fragmentation is important in the correlation of nit with (L^y) since the momentum of 
the b hadron determines the relativistic boost in Eq. 0] Thus, uncertainties in the mean, 
(Xb), and in the shape of the Xb distribution could contribute systematic uncertainties 
to rrit. Since (L^y) is directly proportional to the average transverse momentum of the b 
hadrons, we determine the uncertainty associated with (Xb) by simply varying (L^y) by 
a fractional amount equal to the fractional uncertainty in (Xb). Variations of the shape 
of the Xb distribution could affect (L^y), even for a fixed {Xb), through the processes of 
event selection and b tagging. Nevertheless, since OPAL data does not strongly distinguish 
between a variety of theoretical models j2l| for b fragmentation in Z decays, we do not 
associate any uncertainty with the shape of the fragmentation function. 

A potential source of uncertainty associated with our treatment of b fragmentation stems 
from the fact that the color environment in tt events is different from that of Z events. 
However, recent calculations j^] indicate a value of (Xb) in top quark decays that is within 
one standard deviation of the value in Eq. IHl We have therefore assumed that the results of 
previous measurements can be applied to tt events. 

Taking the preceding discussion into account, we estimate a systematic uncertainty in 
rrit associated with b fragmentation by varying (L^y) by ±0.57%, corresponding to the total 
uncertainty in Eq. IHl 

The fragmentation process results in the production of all b hadron species. These appear 
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TABLE I: Fractions of b hadron species in Z ^ bb decays and their lifetimes. 



^T^f^m PC! 




T ,1 TP'I'i Tn p 






39.7 ± 1.0 


1.536 ±0.014 


1.0 




39.7 ± 1.0 


1.671 ±0.018 


1.0 


Bs, Bs 


10.7 ± 1.1 


1.461 ±0.057 


1.0 


b mesons 


90.1 ± 1.7 


1.603 ±0.022 


1.0 


b baryons 


9.9 ± 1.7 


1.208 ±0.051 


0.8 ±0.05 



in specific proportions, as measured by CDF and the LEP and SLD experiments 26|, and 
with specific hfetimes As for h fragmentation, uncertainties in the h hadron hfetimes 
are important since {L^y) depends on them directly. To take this into account in our studies. 



we use the current uncertainty on the inclusive h hadron lifetime 

1.574 ±0.008 X 10"^^ 



id 



n 



sec 



(7) 



This corresponds to a 0.51% variation in {L^y) if all types of h hadrons are tracked and 
vertexed with equal efficiency. If this is not the case, then the variation may need to be 
amended and an uncertainty must be introduced in association with the modeling of h hadron 
decays. 

The current values and uncertainties for the proportions and lifetimes of the various h 
hadron species are summarized in Table HI For Monte Carlo tt events with full simulation of 

n 

the CDF detector and h tagging algorithm |1^, the h meson species are tagged at a uniform 
rate, while the h baryons tagging efficiency is 80% of the meson efficiency [22|. This means 
that the effective mean lifetime of the h hadrons entering into our {L^y) estimate is higher 
than the average h hadron lifetime in Eq. [7| by an amount Ar = 0.007 ± 0.002 x 10~^^. 
This represents a 0.40 ± 0.1% shift where the error takes into account the uncertainties in 
tagging and in the h baryon fraction, but does not include an uncertainty associated with the 
modeling of b baryon decays, for which there exists limited data at this time. For this reason 
we inflate our estimate by a factor of two, resulting in a 0.2% uncertainty in {L^y) resulting 
from the lower h tagging efficiency for h baryons. It is our expectation that Tevatron Run 
II studies of h baryons will permit a more rigorous treatment of this uncertainty. 
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C. Uncertainties associated with Monte Carlo modeling of tt events. 

There are a variety of uncertainties associated with the PYTHIA Monte Carlo modeling of 
ti events that pertain to the production of jets. To ascertain the corresponding uncertainties 
in nit, we follow the direct advice of the authors [3| or the current conventions of the 
CDF top quark mass group to vary control parameters involved in ti event generation. 
The resulting datasets are then treated in the same manner as the default sample, and 
the resulting (L^y) values are used to extract values of rrit by means of the default mass 
estimator. The difference between the value obtained for each variation and that obtained 
with the default sample is then taken to be the corresponding systematic uncertainty. 

The aspects of ti event generation that could affect jet production rates, in order of 
appearance in the Monte Carlo generation process are: 

1. Spin correlations. 

2. Initial state radiation. 

3. Final state radiation. 

4. Interference of initial and final state radiation. 

5. Multiple interactions. 

Many of these are constrained by experimental results. Final state radiation inW —>■ qq' 
decay, for instance, is equivalent to that in Z ^ qq decay which is very well understood 
from LEP I Q- 

Spin correlations in PYTHIA are present for the b and the W decay products, but not 
for correlations between the two top quarks. In real experiments, there are also multiple 
interactions in beam-beam crossings. These phenomena are expected to have small effect 
on rrit results for this method and were not included in our studies. 

Lowest order Feynman diagrams for initial and final state radiation from top quarks 
are shown in Fig. IHl To estimate the uncertainties associated with final state radiation we 
follow guidance provided by the authors of PYTHIA to the CDF top group j3l in regard to 
the variation of the scale, Aqcd, and the degree of interference with initial state radiation. 
For initial state radiation, this guidance is again applied to our LHC study, while for the 
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FIG. 6: Feynman diagrams for lowest order initial and final state radiation in ti production. 
Interference between these processes is also taken into account. Final state radiation is also possible 
from the b quark and from the quarks appearing in hadronic decay of the W. 

Tevatron, we use variations currently used by the CDF experiment that were determined by 
their studies of Drell-Yan di-muon events j30|. 

With regard to final state radiation, it is thought that PYTHIA may underestimate the 
gluon emission rate relative to the ti + jet matrix elements at LHC energies Q] and so 
one might suspect, a priori, that the prescription used for the Tevatron could be inadequate 
for estimating LHC uncertainties. However, the event selection criteria that we use for the 
LHC, as presented below, is specifically chosen to suppress events with significant radiation. 
We therefore used the same parameter variations for final state radiation for both cases. 

For all of our studies we use the CTEQ5 IsiJ parton density functions (PDFs) to create 
our default Monte Carlo event samples. To estimate the uncertainties associated with the 



PDFs, we generate ti events using two variations of the MRS-98 [3^ PDFs for which Aqcd, 
has values of 300 GeV and 229 GeV. The difference in final (L^y) values between the default 
CTEQ5 sample and the MRS-98 sample with A = 300 GeV is added in quadrature with the 
difference in the two MRS-98 samples resulting from the variation of Aqcd- This procedure, 
while ad hoc in appearance, follows historical conventions at hadron colliders. Recently, 

lowever, a more sophisticated alternative has been suggested by the CTEQ collaboration 

3^ which we will apply in section HVl 

To estimate an uncertainty on the modeling of the top quark spectrum at the Tevatron, 
we took the tree-level spectrum obtained from PYTHIA and re-weighted the events to match 
the NNLO spectrum j2J] at 1.96 TeV, normalized to equal area. We found the resultant 
shift in the extracted top quark mass to be 0.9 GeV/c^ for a generated mass of rrit = 178 
GeV/c^. We draw several conclusions from this. First, for the purpose of our studies, the 
difference is small enough to justify our use of PYTHIA without inclusion of higher order ma- 
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trix elements. Second, the difference between the tree-level and NNLO spectra is expected 
to be significantly larger than the uncertainty on the NNLO calculation alone. Hence, in 
an actual measurement of m^, one would avoid a significant uncertainty by using a Monte 
Carlo for which the top quark px spectrum is consistent with the NNLO prediction. Fur- 
thermore, at the LHC the ti production rates are large enough to allow one to define an 
unbiased sample from which one could do a direct top quark px spectrum measurement. 
The Monte Carlo could be compared to measurement and tuned accordingly to allow the 
associated uncertainty in rrit to be reduced to a negligible level. We therefore do not assign 
an uncertainty associated with the top quark pt spectrum in our studies. 



III. ESTIMATES BASED UPON CURRENTLY AVAILABLE INFORMATION. 

As seen in Section Inl many of the systematic uncertainties associated with this method 
are largely determined by the precision of experimentally determined quantities, such as the 
average b hadron lifetime. In this section we present estimates for the uncertainty on rrit 
for the Tevatron and the LHC using only information that is available at the present time. 
This has the virtue of yielding very sound estimates of systematic uncertainties. It does not, 
however, give a true picture of the precision that might be achieved with this method in 
the future, when much of the information that the method relies upon will be more refined 
that it is at present. This is particularly true in regard to the application of this method 
at the LHC where direct measurements with LHC tt data will allow several of the larger 
systematic uncertainties to be reduced. This is discussed in Section HVl below. 

As mentioned in section |ll the differences between the two environments lead to different 
event selection criteria. We focus on the lepton + jets event signature for the Tevatron, 
and dilepton event signature for the LHC, (see Fig. [7j). In both cases, the term lepton 
actually refers to either a muon or electron from W decay and not a tau lepton, which is 
most often manifested experimentally as a narrow hadronic jet. Events in which the W 
decays leptonically tend to have smaller backgrounds, and hence dilepton event samples, in 
which both W^s decay leptonically, are the most pure. The all-hadronic ti events are those 
for which both W^s decay to quarks. The dilepton, lepton+jets, and all-hadronic ti events 
represent roughly 5%, 30%, and 45% of all ti events, respectively. 

At the Tevatron, where ti yields are much smaller than at the LHC, the lepton+jets 
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FIG. 7: The final state signatures for tt events are determined by the W, which can decay either 
to quarks or to leptons. To determine mt at the Tevatron, we studied the lepton + jets signature 
for ti events in which one W decayed leptonicahy and the other hadronicahy. For the LHC we 
studied the dilepton signature for ti in which the W^s decay to an electron or muon and their 
respective neutrino. 

category provides the most significance. Our Tevatron study, therefore, is focused on this 
signature. We do, however, estimate the extent to which the uncertainty in can be 
reduced by inclusion of dilepton events. For Tevatron event selection criteria, backgrounds 
contribute an uncertainty to rrit as a result of the uncertainties associated with both the 
normalization and shape of the L^y distribution for each of the various background types. 

At the LHC, it is more important to minimize systematic uncertainties. Thus, we focus 
on the dilepton signature to suppress backgrounds and reduce to minimize the occurrence 
of QCD radiation from final state partons. For the same purpose, we normalize our event 
counts to be consistent with the requirements that both b jets are tagged, and that there 
are no other jets in the event. As mentioned earlier, the "jet veto" also helps to minimize 
the probability that there is significant QCD radiation in the event. 

A. Estimates for the Fermilab Tevatron 

Run II of the Tevatron is currently in progress. At the time of this writing, data from an 
integrated luminosity ~0.5 fb~^ have been accumulated by each of the two collider experi- 
ments, CDF and DO. By the end of Run II it is expected that this will reach a value of 4.5 
to 8.5 fb-^ 
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As mentioned above, the use of the lepton+jets signature for tt events is a compromise 
between sample purity and tt event counts. Uncertainties associated with the modehng and 
normahzations of backgrounds contribute a Tevatron-specific uncertainty in rrit. We thus 
start with a discussion of our treatment of these backgrounds before summarizing all of our 
projected systematic and statistical uncertainties. 

1. Tevatron event selection. 

For our modeling of ti events at the Tevatron, the tracks used to determine whether a 
given 6 or c hadron decay could produce an observable displaced vertex, as discussed in 
section Hi B ![ have impact parameter resolution comparable to that of the CDF experiment 
jigj ]. Thus, in Eq. we use parameters: a = 36 /im and /5 = 25 /im-GeV/c. Events are 
selected by requiring one lepton and one neutrino with transverse momentum of at least 20 
GeV/c. There must be a total of 3 or more well-separated partons with transverse momenta 
of at least 15 GeV/c. 

2. Tevatron backgrounds. 

We considered various processes that can have similar characteristics to ti events, and 
thus represent background constituents of Tevatron data sets. The relative proportion of any 
particular background type depends strongly on the event selection. We employ selection 
criteria used by the CDF collaboration for the lepton + jets signature and therefore 
also assume the same types and relative proportions of backgrounds, with the exception of 
WW, WZ, ZZ and Z ^ tt events which together represent only 2.3% of all backgrounds 
and so were ignored for simplicity. We thus consider the following final states: Wbb, Wcc, 
Wc, bb + jets, and mistags in VT+jet events. Background sources are listed in Table HTl 

The contributions from single top quark production were also investigated. The t-channel 
production of a single top quark involves a rather different px spectrum than that found in 
ti events and could therefore be treated as a "background" . However, it makes a negligible 
contribution after our event selection criteria are applied. The s-channel production of a 
single top quark has essentially the same dependence of {L^y) on mt as that of ti events and 
thus represents a potential signal source, but turns out to be negligible in comparison to the 



23 




FIG. 8: The Lxy distribution for ti with mt = 175 GeV/c^, and the L^y distributions for 
various backgrounds in lepton + jets events at the Tevatron. Default Monte Carlo parameters 
were used in addition to the parametric detector resolution model discussed in the text. The signal 
to background ratio is 3 to 1. 



TABLE II: Various sources and proportions of backgrounds used in our studies, following ref. 1 



Source 


Fraction 


Wbb 


20.5% 


Wcc 


7.7% 


Wc 


8.9% 


bb + jets 


30.4% 


PT+jets (mistag) 


32.5% 



contribution from tt. For simplicity, we do not include either process in our studies. 

Backgrounds were .„.ulated using ALPGEN matrix elements Q, (deluding up to 2 add. 
tional radiated partons), together with the HERWIG parton shower model [15|]. Each sample 
was processed with the same selection criteria as that applied to the tt signal. For back- 
grounds involving charm, the c hadron daughter tracks were used for the tagging simulation. 
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and the kinematic cuts previously applied to b jets were applied to the c jets. 

The signal for a given rrit was combined with the total background in a 3-to-l ratio. The 
L^y distributions for tt, {rrit = 175 GeV/c^), and backgrounds is seen in Fig.|Hl The resulting 
values of (L^y) for the various rrit choices were then used to define a mass estimator for the 
Tevatron study. 

3. Tevatron Results. 

We determine the statistical uncertainty as described in section |n] for two different in- 
tegrated luminosities at the Tevatron: the current value of 500 pb~^ and the maximum 
projected Run II value of 8.5 fb~^. For each case, we estimate the corresponding numbers of 
top quarks and background events that we would expect to have by scaling results recently 
reported by CDF For top quark masses in the range 130 < rrit < 230 GeV/c^, we 

create ensembles of large numbers of Monte Carlo samples as described in section m and 
determine {Lr^y) for each. The mass estimator is then used to convert each value of {Lxy) to 
a value of rrit. The distribution of "observed" values of rrit for a given ensemble has average 
value consistent with the mass value input to the Monte Carlo. From the distribution of 
observed rrtt values for a given input mass and choice of integrated luminosity, we determine 
statistical uncertainties by constructing Neymann frequentist confidence intervals. FigureslHl 
andllOlshow the results for integrated luminosities of 500 pb~^ and 8.5 fb~^, respectively. In 
these figures, the dashed lines represent the confidence intervals used to determine statistical 
uncertainties. Thus, for instance, for a top quark mass of rrit = 178 GeV/c^, we expect that 
measurements done at the Tevatron will have statistical uncertainties of 20.8 GeV/c^ and 
5.0 GeV/c^ for integrated luminosities of 500 pb^^ and 8.5 fb~^, respectively. 

The systematic uncertainties, obtained by the methods described in section |nj are the 
same for both cases and are listed in Table ITTTl The statistical uncertainty on the individual 
systematics listed in the table is ~0.2 GeV/c^. 

To estimate the effect of adding tt dilepton events, we scale event counts and tagged h jets 
based upon the results from a recent CDF tt dilepton cross section measurement ^36]. For an 
integrated luminosity of 8.5 fb^^ the statistical error is reduced from 5.0 to 4.4 GeV/c^. For 
the most part, the systematic uncertainties associated with the dilepton signature will be 
similar to those for the lepton+jets signature, except for the uncertainty due to backgrounds 
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TABLE III: Systematic uncertainties for nit = 175 GeV/c^ for the Tevatron. 





T Tn r'pyf p mf / 




fl 7 


r llid;! oLd;Lc i cxU.ld;LlUil 
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x.x 
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U.O 


b liadron lifetime 


1 n 

-L.U 


b fragmentation 


0.9 


b tagging 


0.8 


Backgrounds 


0.3 


Jet Energy Scale 


1.9 


TOTAL 


2.9 



which will be smaller. 

B. Estimates for the CERN LHC 

For our LHC study, anticipated tt production rates allow us to define event selection 
criteria that minimize uncertainties. We select only dilepton events in which the two leptons 
have transverse momenta of at least 35 and 25 GeV/c. We also require two b jets with 
transverse energy above 15 GeV/c and no other clustered energy above 10 GeV. These 
event selection criteria effectively eliminate all backgrounds while minimizing the likelihood 
of QCD radiation. The effect of event selection on the transverse momenta of the top 
quarks is seen in Fig ^2 For an integrated luminosity of 10 fb~^, these requirements yield 
an expectation of ~ 13,000 events, out of a ~ 8,000,000 produced tt events. It is estimated 
that these events have a ratio of ti signal to backgrounds in excess of 30-to-l 0- We 
therefore do not include an uncertainty associated with backgrounds in our LHC study. 

For our modeling of tt events at the LHC, we use track impact parameter resolution 
parameters comparable to those of the CMS experiment js^: a = 15 /xm and /3 = 75 
/xm - GeV/c. 

We estimate a statistical uncertainty of 0.9 GeV/c^ for the LHC in 10 fb^^, as seen 
in Fig El The systematic uncertainties are listed in Table IIVI As for the Tevatron, the 
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FIG. 9: The correlation of the input, ("true"), value of mt with the measured, ("observed"), 
value of rrit (solid) is plotted for ensembles of Monte Carlo lepton + jets events for the Tevatron. 
The dashed lines are the one standard deviation contours (statistical only) for 500 pb~^ integrated 
luminosity. Default Monte Carlo parameters were used in addition to the parametric detector 
resolution model discussed in the text. The signal to background ratio is 3 to 1. 

statistical uncertainty on the individual entries in the table is ~0.2 GeV/c^. 

We did not include an uncertainty associated with b tagging efficiency. It is our expecta- 
tion that the large ti samples available at the LHC will allow the b tagging efficiency to be 
measured as a function of the summed pt of the tracks, or the of the b jet, with small 
uncertainty. To estimate the effect of the jet energy scale uncertainty, we varied iet energies 
by ±3% in accordance with the CMS hadron calorimeter technical design report 



IV. FUTURE IMPROVEMENTS. 



In the preceding sections we presented our estimates for the uncertainties on using 
the mean decay length of b hadrons in ti events at the Tevatron and LHC. We used only 
currently available information for modeling of ti events. For the most optimistic Tevatron 
Run II integrated luminosity scenario, and using events with the lepton + jets signature. 
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FIG. 10: The correlation of the input, ("true"), value of nit with the measured, ("observed"), 
value of mt (solid) is plotted for ensembles of Monte Carlo lepton + jets events for the Tevatron. 
The dashed lines are the one standard deviation contours (statistical only) for 8.5 fb~^ integrated 
luminosity. Default Monte Carlo parameters were used in addition to the parametric detector 
resolution model discussed in the text. The signal to background ratio is 3 to 1. 

we obtained statistical and systematic uncertainties of 5.0 and 2.9 GeV/c^, respectively. 
We estimated that the statistical uncertainty can be reduced to 4.4 GeV/c^ by inclusion of 
events conforming to the dilepton signature, without significantly affecting the systematic 
uncertainty. This corresponds to a total uncertainty of ~5 GeV/c^. Since our estimated 
statistical uncertainty is larger than our estimated systematic uncertainty, we do not expect 
this result to be improved without an increase in total integrated luminosity or an increase 
in the event identification or b tagging efficiencies. With regard to the former some small im- 
provement could be obtained by loosening event selection criteria and allowing background 
to rise until an optimum total uncertainty is attained. With regard to b tagging, we have 
used currently achieved efficiencies, (in the range of 25-40%), that could turn out to be con- 
servative. For instance, the CDF experiment has only recently begun to include information 



from the Layer 00 silicon detector 



H 



4l|, which provides track hits at an average radius 



1.5 cm that leads to a substantial improvement in track impact parameter resolution. The 
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TABLE IV: Systematic uncertainties for rrit = 175 GeV/c^ for the LHC study. 
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FIG. 11: The transverse momentum of the tt pair at ^/s = 14 TeV is plotted for all events, and 
for only those events passing the dilepton event selection discussed in the text. The curves are 
normalized to unit area. 

I — I 

DO collaboration is now constructing a similar small radius silicon detector [42] to enhance 
their b tagging. These devices, as well as gradual improvements in the understanding and 
alignment of tracking detectors at the Tevatron will likely result in higher b tagging effi- 
ciencies. Thus, the actual statistical uncertainties on rrit per Tevatron experiment, for any 
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r >4s=14TeV: 
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Top Mass [GeV/c'] 

FIG. 12: An example of the correlation of Lxy with nit for the default Monte Carlo and our 
particular choice of event selection criteria for the LHC. The correlation is well fit by a straight 
line over small mass ranges such as the one pictured here. Each point represents ~ 26,000 b jets 
in ti events. 

particular integrated luminosity, could be lower than the estimates presented above. 

For 10 fb^^ integrated luminosity at the LHC, we projected statistical and systematic 
uncertainties of 0.9 and 2.4 GeV/c^, respectively, corresponding to a total uncertainty of 
2.5 GeV/c^. Such precision is already significant, as it would allow for an important cross- 
check or improvement of other measurements of rrit. There is, however, the possibility for 
reduction of several of the systematic uncertainties associated with our LHC study. 

As seen in Table IIVI the dominant systematic uncertainties at the LHC are those asso- 
ciated with initial state radiation, b fragmentation, and the average b hadron lifetime. The 
uncertainty in the average b hadron lifetime, and lifetimes of the individual b hadron species, 
continue to improve with time and the accumulation of significantly larger bb event samples 
at B factories and the Tevatron. It is therefore reasonable to anticipate that by the time 
this measurement is performed at the LHC, there will be a reduction in the uncertainties 
associated with the b lifetimes. The uncertainties associated with QCD radiation will be 
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FIG. 13: The correlation of the input, ("true"), value of mj with the measured, ("observed"), value 
of rat (sohd) is plotted for ensembles of Monte Carlo dilepton events for the LHC. The dashed lines 
are the one standard deviation contours (statistical only) for 10 fb^^ integrated luminosity. 

significantly improved as a result of direct measurements at the Tevatron and LHC. Event 
samples containing millions of tt events at the LHC will allow stringent constraints to be 
placed on the Monte Carlo modeling of QCD radiation that could reduce the corresponding 
systematic uncertainties to negligible levels. This in turn will allow event selection criteria 
to be loosened considerably to reduce the statistical uncertainty without significant penalty 
in the systematic uncertainty on m^. 

On the other hand, there are other systematic uncertainties that may not improve. For 
example, the average momentum carried by the h hadron in the process of fragmentation 
is an experimentally determined quantity that is currently dominated by the measurements 
with Z ^ bb events at LEP and SLD, which are limited by systematic uncertainties 
These measurements also benefited a great deal from the beam-energy constraint provided 
by collisions of electrons with positions and are thus unlikely to be significantly improved at 
hadron colliders. It may, however, be possible for some overall improvement to be obtained 
by combining existing measurements. 
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With regard to parton distribution functions, it is difficult for us to say how much LHC 
data might improve the estimates we presented in Section llll Bl although some improvement 
seems inevitable. However, as we noted in Section 111 CI a more sophisticated algorithm PDF 
uncertainties has recently been developed. The CTEQ collaboration has established a new 
framework that enables correlations between the experimental measurements used to define 
the PDFs to be taken into account. It provides a pragmatic way to quantify uncertainties via 
an eigenvector approach to the Hessian method j33|. The eigenvectors form an orthonormal 
basis for a 20 dimensional parameter space of the 20 free parameters of the CTEQ6M fit. 
Forty-one PDF sets are generated, corresponding to positive and negative excursions along 
each of the 20 orthogonal axes in the parameter space, and the central CTEQ6M set. These 
can be used to estimate the uncertainty for any physical quantity that depends upon them, 
including the top quark mass, by adding in quadrature the differences between the positive 
and negative sets for each of the eigenvectors. When we apply this procedure to our Tevatron 
analysis, we obtain an uncertainty that is ~15% smaller than our previous result shown in 
Table Uni 

These considerations lead us to conclude that by the time a measurement of the top quark 
mass has been performed at the LHC with 10 fb^^ of integrated luminosity, the result will 
be substantially more precise than the estimate we presented in section IIII Bl For example, 
a factor of two reduction in the uncertainties associated with QCD radiation, together with 
smaller 10-40% reductions in the uncertainties associated with PDFs, b fragmentation, and 
average b hadron lifetime, reduce our total estimated uncertainty to ~ 1.5 GeV/c^ per LHC 
experiment. 

V. CONCLUSION 

We have presented a new method for the measurement of the top quark mass at the 
Tevatron and LHC that is largely uncorrelated with other methods. Initial studies using 
currently available information to model all aspects of ti events yields uncertainties on 
of ~5 GeV/c^ for 8.5 fb~^ of integrated luminosity in Run II of the Tevatron and ~2.5 
GeV/c^ for 10 fb^^ of integrated luminosity at the LHC. Furthermore, with likely improve- 
ments in our understanding of b hadron properties and tt events the uncertainties associated 
with this method at the LHC could be substantially reduced, making this method com- 
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parable in mass resolution to all other methods, for which a total uncertainty of 1.5-2.0 
GeV/c^ per experiment is expected. 
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